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Abstract 

We address the problem of the separation of variables for the Hamilton- 
Jacobi equation within the theoretical scheme of bi-Hamiltonian geome- 
try. We use the properties of a special class of bi-Hamiltonian manifolds, 
called LvN manifolds, to give intrisic tests of separability (and Stackel 
separability) for Hamiltonian systems. The separation variables are nat- 
urally associated with the geometrical structures of the ojN manifold 
itself. We apply these results to bi-Hamiltonian systems of the Gel'fand- 
Zakharevich type and we give explicit procedures to find the separated 
coordinates and the separation relations. 
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1 Introduction 

The technique of additive separation of variables for solving by quadratures the 
Hamilton- Jacobi (HJ) equation is a very important tool in analytical mechanics, 
initiated by Jacobi and others back in the nineteenth century (see, e.g., I^S], 
^). Following these classical works, an n-tuple {Hi,... ,Hn) of functionally 
independent Hamiltonians will be said to be separable in a set of canonical 
coordinates (gi, . . . , g„,pi, . . . if there exist n relations, called separation 
relations, of the form 



(l)i{qi,Pi,Hi, . . . ,Hn) = , i = l,...,n, with det 



^0. (i.i: 



The reason for this definition is that the stationary Hamilton- Jacobi equations 
for the Hamiltonians Hi can be collectively solved by the additively separated 
complete integral 

n 

W{qi, . . . , tti, . . . , a„) = ^ Wi{qi; «!,...,«„), (1.2) 

i=l 

where the Wi are found by quadratures as the solutions of ordinary differential 
equations. 

One of the first systematic results was found by Levi-Civita, who provided, 
in 1904, a test for the separability of a given Hamiltonian in a given system 
of canonical coordinates. Stackel and Eisenhart concentrated on Hamiltonians 
quadratic in the momenta and orthogonal separation variables. In particular, 
Stackel considered the Hamiltonian 
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and showed that H is separable in the coordinates if there exist an invert- 
ible matrix S{q) and a column vector U{q) such that the z-th rows of S and U 
depend only on the coordinate gj, and H is among the solutions (ifi, . . . , Hn) 
of the linear system 



n ^ 

^ Sij{qi)Hj = -pI - Ui{qi) . 



These equations provide the separation relations for the (commuting) Hamilto- 
nians (i^i, . . . , if^). 

With the works of Eisenhart, the theory of separation of variables was in- 
serted in the context of global Riemannian geometry, and this still represents 
an active area of research, where the notions of Killing tensor and Killing web 
play a key role (see, e.g., Q). 

Starting from the study of algebraic-geometric solutions of (stationary re- 
ductions of) soliton equations and the introduction of the concept of algebraic 
completely integrable system 0, |^, separation of variables has received a 
renewed attention (see, e.g., [|^, |l], |3^). This research activity, also 



connected with the theory of quantum integrable systems, deals with Hamilto- 
nian systems admitting a Lax representation with spectral parameter and an 
r-matrix formulation. In this case, the separation relations are provided by the 
spectral curve 

det(/i/-L(A)) = 

associated with the Lax matrix LiX). Indeed, one can often find canonical coor- 
dinates (Ai, . . . , A„, fii, . . . , fin) on the phase space such that every pair (Aj, /ij) 
belongs to the spectral curve. Since the Hamiltonians are defined by the spectral 
curve, they are separable in these coordinates. 

The two classes of separable systems briefly recalled above strongly suggest 
that a "theory of separability" should start from the following data, 

1. A class of symplectic manifolds M; 

2. A class of canonical coordinates on M; 

3. A class of Hamiltonian functions on M, 
and should provide 
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a) separability test(s) to ascertain whether the HJ equations associated with 

the selected Hamiltonians admit a complete integral which is additively 
separated in the chosen coordinates; 

b) algorithms to compute the separation coordinates and to exhibit the sepa- 

ration relations, so that the HJ equations can be explicitly solved. 

In the context of Riemannian geometry, the manifolds are cotangent bundles 
of Riemannian manifolds, the coordinates are (fibered) orthogonal coordinates, 
and the Hamiltonians are quadratic in the momenta. For Lax systems, roughly 
speaking, the manifolds are suitable coadjoint orbits in loop algebras, the coor- 
dinates are the so-called spectral Darboux coordinates possibly to be found 



using the "Sklyanin magic recipe" p9| , and the separable Hamiltonians are the 
spectral invariants. 

The point of view herewith presented is the following. The class of mani- 
folds we will consider are particular bi-Hamiltonian manifolds, to be termed uN 
manifolds, where one of the two Poisson brackets is nondegenerate and thus de- 
fines a symplectic form uj and, together with the other one, a recursion operator 
N. The class of coordinates, called Darboux- Nijenhuis (DN) coordinates, are 
canonical with respect to uj and diagonalize N . 

The first result is that an n-tuple {Hi, . . . , Hn) of Hamiltonians on M (where 
n = I dim M) is separable in DN coordinates if and only if they are in involution 
with respect to both Poisson brackets. This condition is clearly intrinsic, i.e., 
it can be checked in any coordinate system. A second result of the present 
paper is that examples of separable systems on ujN manifolds are provided by 
suitable reductions of bi-Hamiltonian hierarchies, called Gel'fand-Zakharevich 
systems. They are bi-Hamiltonian systems defined on a bi-Hamiltonian manifold 
(M, {■,■},{■,■}') by the coefficients of the Casimir functions of the Poisson pencil 
{■,-}x '■= {' r}' — - r} ■ Such coefficients are in involution with respect to both 
Poisson brackets, and are supposed to be enough to define integrable systems 
on the symplectic leaves of {■, ■}. If there exists a foliation of M, transversal to 
these symplectic leaves and compatible with the Poisson pencil (in a suitable 
sense), then every symplectic leaf of {■, ■} becomes an ujN manifold, and the 
(restrictions of the) GZ systems naturally fall in the class of systems which are 
separable in DN coordinates. For this reason, we can say that the Poisson pencil 
separates its Casimirs. 

The third result concerns the Stackel separability. With a slight extension 
of the classical notion, we say that {Hi, . . . ,Hn) are Stackel separable if the 
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separation relations are affine in the Hi. 

n 

'^Sij{qi,Pi)Hj -Ui{qi,pi) = , i = l,...,n. (1.3) 

i=i 

In this case, the collection {Hi,... ,Hn) is called a Stackel basis. We give 
an intrinsic test for the Stackel separability in DN coordinates, which has a 
straightforward application to GZ systems. This goes as follows. We notice 
that if {Hi, . . . , Hn) are in involution with respect to both Poisson brackets (and 
therefore separable in DN coordinates), then there exists a matrix F (depending 
on the choice of the Hj) such that 

n 

N*dHi = J2 FijdHj . 
We prove that {Hi, . . . , Hn) is a Stackel basis if and only if 

n 

N*dF,^ = Y,F^kdFuj . 

k=l 

The geometric theory of separability we present in this paper may be, in 
our opinion, regarded as an effective bridge between the "classical" and the 
"modern" aspects of the theory of separability. More evidence of this claim will 



be given in where we will also show how to frame Eisenhart's theory within 
our approach, and discuss the problem of associating a Lax representation to 
GZ systems. 

This paper is organized as follows. The first part (Section 2 to 5) is devoted 
to the geometry of separability on uN manifolds. In Section |^ we will introduce 
the notion of uN manifold and we will study the DN coordinates. Section ^ 
contains the main results about separability on ujN manifolds, whereas in Sec- 
tion H the Stackel separability is considered. In Section ^ we will come back to 
DN coordinates, pointing out some algorithms for their explicit computation. 

In the second part of the paper we will turn our attention to GZ systems. 
Section ^ deals with the particul where there is only one Casimir of the 

Poisson pencil (i.e., one bi-Hamiltonian hierarchy), and contains the example 
of the 3-particle open Toda lattice. This section is intended for an introduction 
to Section |^, where the general case is treated. We will give conditions under 
which a bi-Hamiltonian manifold is foliated in uN manifolds, and we will show 
that the GZ systems are separable in DN coordinates. Subsection |7jB| is devoted 
to the Stackel separability of such systems. In Section we will show an efficient 
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way to determine, in the Stackel separable case, the separation relations for GZ 
systems. Finally, we present an example in the loop algebra of s[(3). 

Acknowledgements. The results presented in this paper are a first account of 
a long-standing collaboration with Franco Magri, which we gratefully acknowl- 
edge. We wish to thank also Sergio Benenti, Boris Dubrovin, and John Harnad 
for useful discussions. This work has been partially supported by INdAM- 
GNFM and the Italian M.U.R.S.T. under the research project Geometry of 
Integrable Systems. 



2 ujN manifolds 

In this section we describe the manifolds where our (separable) systems will be 
defined. They are called uN manifolds, since they are Poisson-Nijenhuis (PN) 
manifolds |29| such that the first Poisson structure is nondegenerate. 



and therefore defines a symplectic form. In turn, PN manifolds are particular 
instances of bi-Hamiltonian manifolds, i.e., smooth (or complex) manifolds M 
endowed with a pair of of compatible Poisson brackets, {■, ■} and {■,•}'. This 
means that every linear combination of them is still a Poisson bracket. 

Definition 2.1 An ujN manifold is a bi-Hamiltonian manifold (M, {■,■},{■,■}') 

in which one of the Poisson brackets (say, {■,■}) is nondegenerate. 

Therefore, M is endowed with a symplectic form u defined by 

{f,g}=ujiXj,X,) , (2.1) 

where Xf is the Hamiltonian vector field associated with / by means of {■, ■}. 
In terms of the Poisson tensor P corresponding to {■,■}, viewed as a section of 
Hom(T*M, TM), this simply means that P is invertible and u is its inverse. 
Using also the Poisson tensor P' associated with {■,■}', one can construct the 
tensor field N := P'P^^, of type (1, 1), to be termed recursion operator of the 
ujN manifold M. 

Proposition 2.2 The Nijenhuis torsion of N, 

T{N){X, Y) := [NX, NY] - N{[NX, Y] + [X, NY] - N[X, Y]) , (2.2) 
vanishes as a consequence of the compatibility between P and P' . 
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A proof of this well known fact can be found in . 

There are two main sources of examples of ujN manifold. The first one comes 
from classical mechanics. Let Q be an n-dimensional manifold endowed with 
a (1, 1) tensor field L with vanishing Nijenhuis torsion, and let us consider its 
cotangent bundle T*Q with the canonical Poisson bracket {■,■}. As shown in 



23| , the vanishing of the Nijenhuis torsion of L entails that one can use it to 



define a second Poisson bracket {■,■}' on T*Q as 

^ - ^ j Pfe > 

where {pi, qi) are fibered coordinates. This Poisson bracket is compatible with 
{■,■}, so that the phase space T*Q becomes an ljN manifold, whose recursion 
operator is the complete lifting of L. 

The second class of examples of uN manifolds can be obtained by reduc- 
tion from a bi-Hamiltonian manifold (M, P, P') where both Poisson tensors are 
degenerate (see, e.g., |ll4[). 

This happens, in particular, in the following situation. Suppose that P has 
constant corank k, that dimM = 2n + k, and that one can find a /c-dimensional 
foliation Z of M with the properties: 

1. The foliation Z is transversal to the symplectic foliation of P; 

2. The functions which are constant along Z form a Poisson subalgebra of 
(C°°(M), {-, ■}) and of (C~(M), {-, ■}'), i.e., if / and g are constant along 
Z, then the same is true for {/,(?} and {/,(?}'. 

Then any symplectic leaf S of {■,■} inherits a bi-Hamiltonian structure from 
M. Moreover, the reduction of the first Poisson structure coincides with the 
symplectic form of S, so that 5* is an ujN manifold. Such a procedure is one of 
the main topics of the paper, and will be fully discussed in Section |^, where we 
will also show that bi-Hamiltonian systems on M give rise to separable systems 
on 5*. The corresponding variables of separation are going to be introduced in 
the next subsection. 

2.1 Darboux— Nijenhuis coordinates 

In this subsection we will describe a class of canonical coordinates on uN mani- 
folds, called Darboux-Nijenhuis coordinates. They will play the important role 
of variables of separation for (suitable) systems on uN manifolds. 
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Definition 2.3 A set of local coordinates {xi,yi) on an ujN manifold is called 
a set of Darboux-Nijenhuis (DN) coordinates if they are canonical with respect 
to the symplectic form u, 



n 

UJ 

1=1 



dyi A dxi , 

1=1 

and put the recursion operator N in diagonal form, 

^ = 5Z (^a^ ® c^Xi + — ® dyij . (2.3) 
This means that the only nonzero Poisson brackets are 



The assumption, contained in ( p.3|) , that the eigenvalues Aj of N are (at least) 
double is not restrictive, since its eigenspaces have even dimension, equal to the 
dimension of the kernel of P' — XiP. For the uN manifold T*Q described in 
the previous section, it is easy to check that the eigenvalues of L (if they are 
independent) and their conjugate momenta are DN coordinates. In order to 
ensure the existence of DN coordinates on more general uN manifolds, we give 
the following 

Definition 2.4 A 2n- dimensional uN manifold M is said to be semisimple if 

its recursion operator N has, at every point, n distinct eigenvalues Ai, . . . , A„. 
It is called regular if the eigenvalues of N are functionally independent on M . 



It can be shown |T8|, |2^, ^ that every point of a semisimple uoN manifold has a 
neighborhood where DN coordinates can be found, and that, if the ujN manifold 
M is regular, one half of these coordinates are "canonically" provided by the 
recursion operator. Indeed, as a consequence of the vanishing of the Nijenhuis 
torsion of A^, the eigenvalues Aj always satisfy 



N*d\i = X.dXi 



I 1 



where A^* is the adjoint of N, and one has 

Proposition 2.5 In a neighborhood of a point of a regular ujN manifold where 
the eigenvalues of N are distinct it is possible to find by quadratures n functions 
/ii, . . . , /i„ that, along with the eigenvalues Ai, . . . , A„, are DN coordinates. 
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Such coordinates will be called a set of special Darhoux-Nijenhuis (sDN) coor- 
dinates. They will often be used in the sequel, because the Aj are simply the 



roots of the minimal polynomial of A^. Proposition means also that every 
regular ujN manifold is locally equal to the "lifted" ujN manifold T*Q we have 
seen in Section |[ 

The distinguishing property of the pairs of DN coordinates (xi, ^j), and, a for- 
tiori, of the "special" pairs (Aj,/Xj), is that their differentials span an eigenspace 
of iV*, that is, satisfy the equations 

N*dxi = Xidxi , N*dyi = Xidyi , i = 1,. . . ,n . (2.4) 

This leads us to the following 

Definition 2.6 A function f on an uoN manifold is said to he a Stackel function 

(relative to the eigenvalue Aj of N) if 

N*df = Xidf . (2.5) 

The following property of Stackel functions, which also explains their name, will 
be used many times in the rest of the paper. 

Proposition 2.7 Let M be a semisimple ujN manifold. A function f on M is 
a Stackel function relative to Aj if and only if, in any (some) system {xi, . . . ,yn) 
of DN coordinates, f depends only on Xi and yi. 

Proof. It is obvious that if / = f{xi,yi) then N*df = Xidf. Conversely, if 



( |2.5|) holds, then df belongs to the Aj-eigenspace of A^*, so that df is a linear 
combination of dxi and dyi and therefore / depends only on Xj and y^. 

QED 



3 Separability on ujN manifolds 

In Section 2 we have introduced a class of (symplectic) manifolds and we have 
selected a class of (canonical) coordinates on such manifolds. Now we are going 
to characterize, from a geometric point of view, those integrable Hamiltonian 
systems on ujN manifolds which are separable in DN coordinates. In the next 
section we will consider the same problem for Stackel separability. 
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We recall that an n-tuple {Hi, . . . , Hn) of functionally independent Hamil- 
tonians on an uN manifold M is said to be separable in the DN coordinates 
(xi, . . . , Xn, 1/1, ... , Un) if there exist relations of the form 

d(j)i 



{xi,yi,Hi,...,Hn) = , i = l,...,n, with det 



^0. (3.1) 



It can be easily shown (e.g., via the Hamilton- Jacobi method) that this entails 
the involutivity of the Hi. Obviously enough, the separability property is not 
peculiar of the specific choice of the functions H^. If Ki = Ki[Hi, . . . , if„) are 
functions of the Hi, they are also separable according to ( |3.1|) . So we see that 
the property ( |3.1] ) concerns the geometrical features of an integrable system, 
i.e., is to be regarded as a property of the Lagrangian distribution defined by 
the mutually commuting functions Hi. Thus one can say that the Hi define a 
separable foliation of M. According to the following theorem, that will be proved 
during this section, the separability property can be formulated in terms of the 
geometric objects u and A^, or {■, ■} and {■,■}', of the uN manifold M. 

Theorem 3.1 Let M be a semisimple uoN manifold and let {Hi,... ,Hn) be 
a set of n functionally independent Hamiltonians on M . Then the following 
statements are equivalent: 

a) The foliation defined by {Hi, . . . , Hn) is separable in DN coordinates (and 
therefore Lagrangian with respect to uj ); 

b) The distribution tangent to the foliation defined by {Hi, . . . ,Hn) is La- 
grangian with respect to uj and invariant with respect to N; 

c) The functions {Hi,... ,Hn) are in bi-involution, i.e., {Hi, Hj} = and 
{Hi,HjY = for all 1, 3. 

We will often refer to property c) by saying that the foliation defined by the Hi is 
bi- Lagrangian. This is a fundamental property in our approach to separability, 
and will be exploited especially in Sections |^ and ^ Incidentally, we notice that 
bi-Lagrangian foliations play an important role in the study of special Kahler 



manifolds 21 



Throughout the rest of the section M will be a semisimple ujN manifold, 
(Ai, . . . , A„) the eigenvalues of the recursion operator N, and {xi,yi) DN coor- 
dinates on M. We begin with showing that the invariance with respect to is 
a necessary condition for separability. 

Proposition 3.2 Let {Hi, . . . , Hn) be functions on M that are separable in DN 
coordinates. Then the subspace spanned by {dHi, . . . ,dHn) is invariant with 
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respect to N* . More precisely, there exists a (simple) matrix F with eigenvalues 
(Ai, . . . , An) such that 

n 

N*dHi = FijdHj , i = 1, . . . , n . (3.2) 

Consequently, the Lagrangian distribution defined by {Hi,...,Hn), which is 
spanned by the Hamiltonian vector fields Xh^, is invariant with respect to N . 

Proof. Differentiate the relations ( |3.1|) , 

d(t)i, , d(pi ^ d(j)i 

— + — % + X:^rfi7, = 0, (3.3) 

then apply A^* to obtain 

^A.....|iAA + |:^iVW,^0. (3.4) 



It follows that 



j^^dH, \dxi dvi J j^.dHj 



(3.5) 



that is, in matrix form 



JN*dH = AJdH , (3.6) 

where J^j = §f-, dH = {dHi, . . . , dHnf , N*dH = {N*dHi, . . . , N*dHnf, and 
A = diag(Ai, . . . , A„). Therefore ( p.2| ) is satisfied with F = J^^AJ, and the 
eigenvalues of F are (Ai, . . . , A^). The final assertion easily follows. 

QED 

The matrix F will be called the control matrix, with respect to the basis 
{Hi, . . . , Hn), of the separable foliation. 

Proposition 3.3 // {Hi, . . . , if„) define a distribution which is invariant with 
respect to N , that is, 

n 

N*dHi = FijdHj , i = l,...,n , (3.7) 

and the eigenvalues of F are distinct, then the Hi are separable in DN coordi- 
nates. 
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Proof. Since the eigenvalues of F are distinct, they are the eigenvalues (Ai, . . . , 
of A^, so that there exists a matrix 5* such that F = S^^AS. With S we define 
the 1-forms 9i := J2^=i SijdHj, for i = 1, . . . ,n. They are eigenvectors of N*, 
since 

n n n 

N*di = J2 S^,N*dH.j = SijFjkdHk = J2 >^iSikdHk = XA . (3.8) 

j=l j,k=l k=l 

Then there exist functions Li and Mj such that 9i = Lidxi + Midyi, that is, 

n 

^ SijdHj - Lidx, - Midyi = . (3.9) 

3=1 

This means that dim{dHi, . . . ,dHn,dxi,dyi) < n + 1, so that there exists a 
relation of the form (|3J] ), i.e., the functions {Hi, . . . , Hn) are separable in DN 
coordinates. 

QED 

In order to complete the proof of the equivalence between statements a) and b) 
of Theorem |3.1|, we need the following: 



Lemma 3.4 // {Hi, . . . , Hn) are independent functions in involution with re- 
spect to uj such that ^3.2^ holds, then the eigenvalues of F are distinct. 

Proof. Suppose that N*dHi = Y.]=iFijdHj, with {Hi,Hj} = for all i,j. 
Since F represents the restriction of A^* to {dHi, . . . ,dHn), it is diagonaliz- 
able. Thus, if Aj would be a double eigenvalue of F, the span {dHi, . . . , dHn) 
would contain the 2-dimensional eigenspace spanned by dxi and dyi. But the 
involutivity of the Hi would entail that {xi,yi} = 0, which is false. 

QED 

Relations ( p.2|) may be called generalized Lenard relations (and the functions 

Hi fulfilling them a Nijenhuis chain, as in [1^), as enlightened by the following 
example. 

Example 3.5 If Hk := ^ti N'' = E"=i then N*dHk = dHk+i for k = 
1,. . . ,n — 1, that is, the Lenard relations P'dH^ = PdH^^i hold. Moreover, 
dHn+i = X]j=i CjdHn+i-j, where A" — X]j=o '^n-jX^ is the minimal polynomial 
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of N. Therefore, condition (|3^ ) is satisfied with 



1 
1 

••• 






1 

Cl 



(3.10) 



Remark 3.6 It is well known that functions Hi satisfying the Lenard rela- 
tions are in involution with respect to both Poisson brackets, and so they pro- 
vide a first instance of correspondence between invariant distributions and bi- 
involutivity, which is at the same time trivial and paradigmatic. 
Indeed, it is trivial from the point of view of the theory of separation of vari- 
ables, since such Hamiltonians are easily seen to depend only on (Ai, . . . , A„) 
if the uN manifold M is regular and semisimple. Then the Hamilton- Jacobi 
equations associated with the Hi are trivially separable in the sDN coordinates 



A„, /ii, . . . , Nevertheless, it is paradigmatic with respect to the 



issues of this paper. Indeed, the Hi (that is, the A,) define a distinguished bi- 
Lagrangian foliation, called principal foliation, which coincides with the canon- 
ical fibration n : T*Q ^ Q of classical phase spaces when T*Q is the uN 
manifold considered in Section ^. However, there are in general bi-Lagrangian 
foliations which are different from the principal one, as we will explicitly see 
in Section |^. We are going to show that such foliations are characterized by 
the invariance with respect to N, so that they give rise to separable systems. 
This means that our theory deals with cases in which the Hamiltonians are not 
simply the traces of the recursion operator. In other words, we will deal with 
cases in which the control matrix F of equation ( p. 21) need not be a compan- 



ion matrix of the form (|3.1CI|) . Accordingly, the separable vector fields we will 
consider are tangent to a bi-Lagrangian foliation, but they are not, in general, 
bi-Hamiltonian. 



Proposition 3.7 Let {Hi, . . . , Hn) be independent functions on M. Then (^■2J 
holds, with a matrix F with distinct eigenvalues, if and only if the functions Hi 
are in hi-involution: 

{Hi, Hj} = {H„ H,}' = , for all 1,3 = I,..., n. (3.11) 
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Proof. We know from Proposition |3.3| that condition (|3.2| ), with a simple 



matrix F, imphes separabihty and therefore involutivity with respect to {■, ■}. 
Moreover, 

{Hi, HjY = {dHi, P'dHj) = {dHi, NPdH^) = {N*dHi, PdHj) 

(3.12) 

= Ylk=l^ik{Hk,Hj} . 

showing that {Hi^Hj}' vanishes as weU. 

Conversely, suppose that {Hi,Hj} = {Hi.Hj}' = for all Then the 
foliation Ti. defined by the Hi is Lagrangian with respect to {■,■}, and 

{N*dHi, PdHj) = {dHi, NPdHj) = {dHi, P'dHj) = . 

Thus, N*dHi belongs, for every i, to the annihilator of {P dHi, . . . , P dHn), 
which is tangent to H, since H is Lagrangian. This shows that (|3.2|) holds, and 



Lemma |3.4| entails that F has distinct eigenvalues. 

QED 

Thus we have proved also the equivalence between b) and c) of Theorem |3TT| . 

Remark 3.8 One could also prove that a function F is separable in DN coordi- 
nates if and only if its Hamiltonian vector field Xp is tangent to a bi-Lagrangian 
foliation Ti. The "if" part of this statement is a simple corollary of Theorem ^.1| . 
Indeed, let Ti. be defined by the functions {Hi, . . . , Hn). Then F is a function 
of the Hi, since the distribution is Lagrangian, and one can find other {n — 1) 
functions K2, ■ . ■ , Kn such that 7i is defined by (F, K2, . . . , Kn). The "only if" 
part of this statement is deeper, and essentially gives rise to the intrinsic picture 
of the Levi-Civita conditions for separability, to be fully discussed in |l^ . 



Summing up, we have proved a criterion for the separability in DN coordi- 
nates, which can be tested without knowing explicitly these coordinates. Indeed, 
the statement ( p.ll| ) can be checked in any coordinate system. An important 



application of this criterion will be given in Section |^, where we will show 
that the bi-Hamiltonian hierarchies on a bi-Hamiltonian manifold admitting a 
transversal distribution with the properties mentioned at the end of Section ^ 
give rise to separable Hamiltonian vector fields on the reduced uN manifolds. 

4 Stackel separability on uuN manifolds 

The separability criteria of the previous section do not give explicit information 
on the form of the separating relations ( |3.1| ). For this reason, in this section 
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we will concentrate on the more stringent notion of Stackel separability. Recall 
that {Hi, . . . , Hn), independent functions on an ujN manifold, were defined to 
be Stackel separable in the DN coordinates (xi, . . . ,?/„) if there exist relations 
of the form (|3.1| ), given by affine equations in the Hj, that is, 

n 

'^Sij{xi,yi)Hj-Ui{xi,yi) = , i = l,...,n, (4.1) 

with S an invertible matrix. In this case, we say that the Hi are a Stackel basis 
of the (separable) foliation. The entries and Ui depend only on Xj and yi, 
i.e., they are Stackel functions according to Proposition pl7| . Usually, S is called 
a Stackel matrix, and U a Stackel vector. Notice that the definition of Stackel 
separability depends on the choice of the Hi defining the Lagrangian distribu- 
tion. Indeed, if {Hi, . . . , Hn) are Stackel-separable, then Ki = Ki{Hi, . . . , Hn), 
for i = 1, . . . ,n, will not, in general, fulfill relations of the form (|4.1|) . A nat- 
ural problem, that will not be discussed in this paper, is to give a geometrical 
characterization of the Lagrangian foliations admitting a set of defining func- 
tions for which Stackel separability holds. Some results in this direction will be 
presented in |I0| . 

Now we will give a necessary and sufficient condition for the Stackel separa- 
bility in DN coordinates of a given ra-tuple {Hi, . . . , Hn) of functions on an ujN 
manifold. We will also show that in this case one can explicitly find the rela- 
tions (|3.1j ) and has useful information to algebraically determine the separation 



variables. 

Suppose {Hi,...,Hn) to be independent functions on a regular semisim- 
ple uN manifold that are Stackel separable in the DN coordinates. Then we 



know from Proposition 3.2 that there exists a control matrix F, with eigenval- 



ues (Ai, . . . , A„), such that N*dH = FdH. Since Proposition entails that 
A^*^^ = AdS and N*dU = AdU, we can show: 

Proposition 4.1 In the above-mentioned hypotheses, the matrix F satisfies 

n 

N*dF = FdF , that is, N*dFij = ^ FikdFkj V i, j = 1, . . . , n . (4.2) 

k=l 

Proof. First we show that F = S^^AS. Indeed, 

SFdH = SN*dH = N* [d{SH) - {dS)H] = N*dU - {N*dS)H 
= AdU- AdS H = ASdH . 
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Then we have 

N*dF = N'^diS-^AS) = N'^i-S-^dS S-^AS + S-^dA S + S-^AdS) 

= -S-^AdS S-^AS + S-^AdA S + S'^A'^dS (4.4) 
= S-^AS{-S-^dS S'^AS + S-^dA S + S'^AdS) = FdF , 
and the proof is complete. 

QED 

Condition ( [4 .21) is also sufficient for the Stackel separability, as shown in the 
following: 

Theorem 4.2 Let {Hi, . . . , Hn) be independent functions, defining a bi-Lagrangian 
foliation on a regular semisimple uN manifold. If the control matrix F fulfills 
1^4^, then: 

1. The left eigenvectors of F, if suitably normalized, form a Stackel matrix. 
More precisely, if S is a matrix such that F = S^-^AS , and such that in 
every row of S there is an entry equal to 1, then S is a Stackel matrix in 
DN coordinates (xi, . . . , ?/„); 

2. The functions {Hi, . . . , Hn) are Stackel separable in DN coordinates. 
Proof. From (|4.2| ) we have that 

N*{-S-^dS S-^AS + ^-^A S + S-^AdS) 

(4.5) 

= S-^AS{-S-^dS S-^AS + S-^dA S + S'^AdS) , 

that is, 

N*{-dS S-^AS + AdS) = A{-dS S-^AS + AdS) , (4.6) 

or {-N*dS+AdS)F = A{-N*dS+AdS). Hence the j-thiow of {-N*dS+ AdS) 
is a left eigenvector of F, relative to Xj. This entails that it is proportional to 
the j-th row of S, i.e., there exists a 1-form aj such that 

ej{-N*dS + AdS) = ajCjS , (4.7) 

where cj is the j-th row vector of the standard basis. Multiplying equation (^4.7| ) 
by e^, where Sjk = ^jSe^. = 1, we obtain aj = 0, so that 

N*dS = AdS . (4.8) 
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In components, this reads N*dSjk = XjdSjk, which imphes (see Proposition! 
that Sjk depends only on xj and yj, i.e., 5* is a Stackel matrix. Finally, the fact 
that U := SH is a Stackel vector follows from 

N*dU = N*{dSH + SdH) = AdS H + SFdH = A{dS H + SdH) = AdU . 

(4.9) 

This completes the proof. 

QED 

The results obtained so far can be summarized in the following statements. An 

n-tuple of functions {Hi, . . . , if„) in involution is separable (in DN coordinates) 
if and only if the span of their differentials is invariant for A^*. Let F be the 
matrix which represents (the restriction of) A^* on such a span. Then equation 
( |4.2| ) represents a test for the Stackel separability of the Hi. Once this test 
is passed, the Stackel matrix is easily constructed as a (suitably normalized) 
matrix that diagonalize F, and the separation procedure can be quite explicitly 
performed. Therefore, in our setting the Hamiltonians provide their Stackel 



matrix as well as the separation relations ([4.1|) . 

We end this section with the following comment on the intrinsic meaning of 
the Stackel separability conditions ( ^l2|) . It is known that, as a consequence 



of the vanishing of the Nijenhuis torsion of A^, the de Rham complex of M 
is endowed with a second derivation rfjv, which is defined to be the unique 
(anti) derivation with respect to the wedge product extending 

dNfiX) = df{NX) = N*df{X) ^^^^^ 

d^e{x, Y) = x{e{NY)) - Y{e{NX)) - e{[x, yu) , 

where / is a function, ^ is a 1-form, X, Y are vector fields on M, and 

[X, YU = [NX, Y] + [X, NY] - N[X, Y] . 

This differential is a cohomology operator (c/^ = 0) and anticommutes with the 
usual exterior derivative d. One notices that the invariance condition ( |3.2|) can 
be equally be written, in matrix notation, as 

dNH = FdH. (4.11) 

Imposing the condition c?^ = on this equation, taking into account the anti- 
commutativity of d and df^f, and translating back dj^f = N*df if / is a function 
on M, one gets 

{N*dF - FdF) AdH = 0. (4.12) 
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So we see that the Stackel separabihty conditions (|4.2|) are nothing but a 
"strong" solution of the equations imposed on the control matrix F by the 
cohomological condition d% = 

5 Special DN coordinates 

In this section we will discuss the problem of explicitly finding sets of special 
DN coordinates on an uoN manifold M. We assume that M be regular and 
complex, so that the eigenvalues (Ai, . . . , A„) of can be used as (half of the) 
coordinates on M. We know that in a neighborhood of a point where the Aj are 
distinct there exist functions (/^i, . . . , forming with the eigenvalues a system 
of DN coordinates, and that the yUj can be computed by quadratures. However, 
they can often be found in an algebraic way, as we will see below. We divide 
our argument in three main points. 

We start remarking that there are simple conditions to be checked the /Xj in 
order to ensure that they form with the Aj a set of DN coordinates. To this 
aim, we observe that the /ij must fulfill two kinds of requirements: 

1. They have to be Stackel functions, that is, they must satisfy N*dfii = 

2. They have to fulfill the canonical commutation relations with respect to 
the first Poisson bracket: {Aj,/ij} = 6ij, = 0. 

In principle, these conditions require the computation of the Aj. We will show 
that this can be avoided, and that a smaller number of equations must be 
checked. The first step is to notice that, once conditions 1 are satisfied, condi- 
tions 2 can be replaced with the n equations 

{Xl + ■■■ + Xn,^^^} = 1 , (5.1) 

which do not require the computation of the Aj, but only of their sum, that is, 
Ci := ^tr N and, consequently, of the Hamiltonian vector field Y := —Pdci = 
Sr=i Indeed, suppose that jjj be a Stackel function, and observe that 

Xi{Xi,fij} = Xi{dXi,Pdfij) = {XidXi, Pdfij) = {N*dXi, Pdfij) 

= {dXi, NPdiij) = {dXi, PN*diij) = {dXi, XjPdfXj) (5.2) 
= Aj{Ai,/ij} , 
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so that {Aj,/ij} = if 2 7^ j. Then equation becomes {Aj,/ii} = 1. In the 
same way one shows that {/Xj, /ij} = 0. Hence, in order to find the /ij coordinate 
we have to look for a Stackel function (relative to Aj) such that (|5.1j ) holds. 

The second point starts from the following idea, which will be extensively 
used in the part of the paper dealing with Gel'fand-Zakharevich systems. Let 
us consider the minimal polynomial 

A(A) = A" - (ciA"-^ + C2A"-2 + ■ ■ ■ c„) (5.3) 

of A^. Using the Newton formulas relating the traces of the powers of and 
the coefficients c, of A (A), one easily verifies that the latter satisfy 

N*dci = dci+i + Cidci , i = 1, . . . ,n — 1 

(5.4) 

N*dCn-l = dCn ■ 

These relations are equivalent to the following equation for the polynomial A (A), 

N*dA{X) = ArfA(A) + A(A)rfci . (5.5) 

Relations of this kind are very interesting for our purposes. For instance, it 
holds: 

Proposition 5.1 Let f{x;X) be a function defined on M, depending on an 
additional parameter A. Suppose that there exists a 1-form af such that 

N*d{f{x; A)) = Xd{f{x; A)) + A(A)a/ . (5.6) 

Then, the function fi defined by fi{x) := f{x;Xi{x)), i.e., the evaluation of 
/(x; A) on X = Xi, is a Stackel function relative to Aj. 

Proof. The differential of fi equals 

df-{x)=df{x;X)\^^^^ + ^^^^\^^^dX, , (5.7) 

where, in the term d {f{x; A)) |;^_;^,, one treats A as a parameter. Applying the 
adjoint of the recursion operator we get 

N*d (fix; A,)) = N*d {fix; A)) |,^,^ + X,^^^^\^^^dX, , (5.8) 
whence the assertion, since A(Aj) = 0. 

QED 



19 



Definition 5.2 We will call a function on M depending on the additional pa- 
rameter X a Stackel function generator if it satisfies ( 15. 6| j with a suitable 1-form 



Lemma 5.3 The space of Stackel function generators is closed under sum and 
multiplication, and is invariant with respect to the action ofY. If f is a Stackel 
function generator and g is a function of one variable, then g o f is a Stackel 
function generator. 

Proof. The only assertion whose proof is not straightforward is the invariance 
with respect to Y. This follows from the fact, already noticed in Example p.5| , 
that y is a bi-Hamiltonian vector field; hence Ly{N*) = Ly{P^^P') = and, 
consequently, Ly(Ai) = 0. 

QED 

It is clear that if fi is a Stackel function relative to Aj for i = 1, . . . ,n, then 

there exists a Stackel function generator /(x; A) such that fi = f{x;Xi), e.g., 
the interpolating polynomial. In terms of the generator, condition (|5.1|) can be 
written as 

r(/(x;A)) = l for A = Ai, i = l,...,n. (5.9) 
For further use, we state and prove the following: 

Proposition 5.4 Let f be a Stackel function generator, and suppose that for 
n > 1 the action ofY on f closes, that is, the relation 

n-l 

i=o 



with Y{aj) = 0, holds. Then equation ^5.{^ can be algebraically solved. 



Proof. There are two cases: in a first instance, suppose that, actually, Y is 
nilpotent, that is, Y^{f) = is satisfied for some n > 1 (whilst y"~^(/) ^ 0). 
Then it is easily seen that y"~^(/)/y"^^(/) is a Stackel function generator ful- 
filling (^. 

On the contrary, if (oq, . . . , On-i) 7^ (0, . . . ,0), then the matrix A represent- 
ing the action of F on $ := (/, Y{f), . . . , Y^~^{f))^ has at least one nonzero 
eigenvalue u, which is a solution of u"- = J2^=o o^ju^ . Let w = {wq, . . . ,Wn-i) 
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be a (left) eigenvector of A relative to u, e.g., the one given by Wn-i 
Wk = u"--''-^ - YTjlo'"^ ak+j+ii^^ for A; = 0, . . . , n - 2. Then 

y{^^ogY,w^Y^{f)^ =1. 

Indeed, X]j=o — ^'^'^ 

Y (m;$) = wA(^ = lyw^ , 

implying ( p. 11 ). 

QED 



= 1 and 
(5.11) 



These arguments reveal a further important aspect of Stackel separability 
within our approach to separation of variables. Indeed, the condition of Stackel 
separability, whose intrinsic form is given by equation ( [4.2| ), entails that the 
matrix of the (suitable normalized) eigenvectors of the control matrix F is a 
Stackel matrix, that is, its columns are Stackel functions of A^*. Since we have 
shown that a way to algebraically find the fii coordinates is to find Stackel 
functions (or Stackel function generators) and to combine them in order to fulfill 
equation (^.9|) , we see that, in the Stackel case, the Hamiltonians themselves 
may algebraically provide the coordinates in which the corresponding flows can 
be separated. 



6 Separability on odd-dimensional bi-Hamiltonian 
manifolds 

This section starts the second (and more applicative) part of the paper, in 
which we will use the results of Sections ^ and ^ to discuss the separability of a 
specific family of integrable systems. They are defined on a class of bi-Hamil- 
tonian manifolds, known in the literature as complete torsionless bi-Hamiltonian 



manifolds of pure Kronecker type (see [19, 34] and the references quoted therein). 



In this section we will consider the simplest case, corresponding to generic 
odd-dimensional bi-Hamiltonian manifolds (while in Section ^ we studied the 
case of regular ujN manifolds, which are generic even-dimensional bi-Hamil- 
tonian manifolds). Their Poisson tensors have maximal rank. The more general 
case will be treated (with detailed proofs) in the next section. 
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Let (M, P, P') be a (2n + l)-dimensional bi-Hamiltonian manifold, and let 
the rank of P be equal to 2n. Suppose that the Poisson pencil P\ := P' — XP 
has a polynomial Casimir function 

n 

h{x) = J2h,x--\ 

This amounts to saying that the functions {Hq, . . . ,Hn), which we assume to 
be functionally independent, form a bi-Hamiltonian hierarchy, starting from a 
Casimir Hq of P and terminating with a Casimir of P', 

PdHo = , P dHi+i = P' dHi , P' dHn = . (6.1) 

In particular, they are in involution with respect to {■, ■} and {■,■}'. If dHo ^ 
at every point of M, then the symplectic foliation of P is simply given by the 
level surfaces of Hq. The restrictions of {Hi, . . . , to a symplectic leaf 5* of P 
form an integrable system (in the Arnold-Liouville sense). The corresponding 
Hamiltonian vector fields are the restrictions to S of Xj := P dHi, where i = 
l,...,n. 

At this point it is natural to wonder whether the bi-Hamiltonian structure 
of M can give information on the separability of the (restrictions of the) Hamil- 
tonians {Hi, . . . , if„). More concretely, one can try to induce an uN structure 
on S in order to apply the separability theorems of Sections |^ and ^. As an- 
ticipated in Section |^, this can be done if there exists a vector field Z which is 
transversal to the symplectic foliation of P and fulfills the following condition: 

C) if F, G are functions on M which are invariant for Z, that is, Z{F) = 
Z{G) = 0, then {F, G} and {F, G}' are also invariant. 

In this case, any symplectic leaf of P inherits a bi-Hamiltonian structure from 
M. Clearly, the first reduced bracket is the one associated with the symplectic 
form of S, so that 5* is an ujN manifold. 

In the following section we will prove that, if Z is normalized in such a way 
that Z{Hq) = 1, condition C) takes the infinitesimal form 

LzP = 0, LzP' = Y NZ , (6.2) 

for a suitable vector field Y . In this case there is a useful form of the reduced 
Poisson brackets {■, and {■, -jl^ on the symplectic leaf S. If / and g are 
functions on S, by definition {/, g)s and {/, g^'g should be computed by taking 
local extensions of / and g which are invariant along Z . But one can avoid the 
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use of invariant functions and consider arbitrary extensions F, G. Then 



{f,g}s = {F,G} 

{f^gYs = {F,Gy + X'{F)Z{G)-X'{G)Z{F), 



where X' := P' dH^ = P dHi and the right-hand sides of the previous equations 
are imphcitly understood to be restricted to S. These equations show that the 
restrictions of (ifi, . . . , Hn) to S are in bi- involution, and then separable in DN 



coordinates because of Theorem |3.1| . We are going to show that they are even 
Stackel separable, by computing their control matrix F and checking that it 
satisfies the condition N*dF = FdF. 

To this purpose, we notice that the Lenard relations ( |6.1|) on M give rise to 
the equations 



N*dH, = dHi+i - Z{Hi)dHi , 
N*dHn = - Z{H^)dHi , 



n — 1 



(6.3) 
(6.4) 



where is the recursion operator of the uN manifold S and " denotes the 
restriction to S. Therefore, the control matrix of {Hi, . . . ,Hn) is given by a 
single Frobenius block: 



-Z{Hi) 1 
-Z{H2) 1 







1 




(6.5) 



So we see that the (restriction to the symplectic leaf S) of the functions q = 
—Z{Hi) are the coefficients of the characteristic polynomial of the matrix F, 
that is, the coefficients of the minimal polynomial of the recursion operator N, 
A(A) = A" — (ciA"~^ + ■ ■ ■ + c„). Recalling that the coefficients of the minimal 
polynomial of N satisfy 



N*dCi = dCi+i + CidCi , N*dCn = CndCn 



(6.6) 



we see that the condition N*dF = FdF for the Stackel separability of the 
Hamiltonians is automatically verified. Hence we have proven 
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Theorem 6.1 The Hamiltonians of a corank-1 torsionless GZ system are Stdckel 
separable in DN coordinates. 

It is worthwhile to notice that the examples previously considered in the litera- 
ture within the theory of quasi-bi-Hamiltonian systems (see, e.g., |32|, ^) fall 
into this class. The link with the classical Stackel-Eisenhart theory of separation 



of variables is discussed in 23 



We remark that the vector field Y appearing in ( |6.2| ) can be chosen to be 
tangent to S. In this case, Y = P d{Z{Hi)) = —Pdci, so that its restriction 
to S is the vector field we used in the previous section to determine the /Xj 
coordinates. (This explains why we made use of the same notation). 

Now we will write the separation equations for the GZ Hamiltonians. The 
Stackel matrix S, being the (normalized) matrix of the left eigenvectors of F, 
is easily seen to be the Vandermonde-like matrix 



S 



A? 



n-l 



A 



n-l 



Ai 1 



A. 



where the Aj are the eigenvalues of N, i.e., the roots of A(A). Therefore, the 
separation relations take the form 



n-2 



+ Hn = Ui{Xi, fXi 



(6.7) 



where (Ai, . . . , A„, /ii, . . . , /i„) are special DN coordinates on 5* and the Ui are 
the entries of the Stackel vector. Such entries can be explicitly computed once 
we have the map sending the DN coordinates to the corresponding point of S, 
as we will check in the example of the 3-particle nonperiodic Toda lattice. 

Another way to arrive at the separation equations is to multiply (|6.3|) by 
A"~' and then to add to (16.41). The result is 



A^* dH{\) = A dH{\) - A{\)dH_ 



1 ) 



meaning that H{X) := Yl^=i Hi\"'~^ is a Stackel function generator according to 
Proposition |0| . Thus, in DN coordinates, H{Xi) = Ui{Xi,fii), which coincides 
with (|6.7| ). We stress that H{X), being a Stackel function generator, can be 
in some cases used to determine the /ij coordinates. Instances of this situation 



are provided by the Toda lattice, as discussed in ||T2[, and by the stationary 
reductions of the KdV hierarchy Here we will present the example of the 
3-particle nonperiodic Toda lattice. 
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Example 6.2 The Hamiltonian of the system is 

3 2 

i^Toda = \Y.P^' + E^^P(^' - ^'^') • (6.8) 
i=l i=l 

As usual (see, e.g., [l^, and [|l6l for the separability), one introduces the 
"Flaschka-Manakov coordinates" (ai, 02, &i, &2, ^s), where 

bi=pi , ai = - exp(g* - g*+^) , 

and consider the manifold M = (C*)^ x C^, or M = R>o^ x R^. We endow it 



with the Poisson pencil P\ = P' — XP given by (see, e.g., and references 
cited therein) 



P, 







aia2 



-aia2 {bi - A) ai (A - 62) ai 







(62 - A) 02 (A - 63) 02 



(A-6i)ai 
(62 - A) ai (A - 62) 02 
(63 - A) a2 



ai 



-ai 



a2 



-a2 



(6.9) 



It has a polynomial Casimir H{X) = HqX^ + HiX + if2, where 
Ho = hi + h2 + h 

Hi = -(6162 + &2&3 + ^3^1 + «1 + 02) 

H2 = bib2h + aih + 0-2^1 ■ 



The Hamiltonian (|6.8|) is related to the coefficients of H{X) by i/xoda = Hi + ^Hq. 

There are two nontrivial flows, given by: 

d d d d d 

Xi = PodHi = ai{bi - b2)- h 02(62-^3)7^ h Oit^ + (02 - 01)7:^ ^27:^- 

oai oa2 obi c'02 ob^, 

d d d 
X2 = PodH2 = ai[a2 + b3{b2 -bi)]— + a2[ai + bi{b3 - b2)]-g^ ^^^^df 

d d 
+ (0163 - a2bi)— + a2bi— . 

db2 obs 

The symplectic leaves of P are the level surfaces of Ho, so that they can be 
parametrized by (ai, 02, 61, 62)- A possible choice for the normalized transversal 
vector fleld is Z = because Z{Ho) = 1 and 



LzP = 



Ly P' = Y KZ 
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with Y = asaf;- Since Y{Ho) = 0, we know that Y = Pd{Z{Hi)) = -Pdci. If 
5* is a symplectic leaf of P, the reduced bi-Hamiltonian structure on S is simply 
obtained by removing the last row and the last column of Py. 









ai 


— ai 











a2 


— ai 











ai 


-a2 









p' 





aia2 
—aihi 
aib2 



—aia2 cti&i 



-a2&2 -di 



-ai&2 

02^2 

ai 




For completeness, we display recursion operator, 



N 



P'sPs' 



hi 


aiipi - 62)702 









62 


-02 








ai/a2 


&i 





-1 


-ai/a2 





&2 



whose minimal polynomial is 

A(A) = A^' + Z{Hi)X + Z{H2) = X^- (61 + 62)A + ai + 6162 • 

The coordinates Ai, A2 are its roots. 

The restrictions of Hi and H2 to the symplectic leaf Hq = c are 

Hi = -Cipi + 62) + + + ^1^2 - cii - ^2 
if2 = C(ai + 6162) - (Ctl + &1&2)(&1 + ^2) + ^2^1 ■ 

We know that H{X) := ^lA + ^2 is a Stackel function generator, and that 
the separation equations are H{Xi) = U{Xi,fii), for i = 1,2. To write them 



explicitly, we need the form of the yUj. They can be found using Proposition q]4 
and the fact that 

Y'{H{X)) = Y{H{X)) . 

This entails that /(A) := logF(if(A)) satisfies F(/(A)) = 1, so that, according 
to the results of Section H, 



fj,i = \ogY{H{Xi)) = log(a26i - 02 Ai 



1,2 



form with the eigenvalues Ai and A2 of a set of (special) DN coordinates. 



26 



Finally, using the expression of (ai, a2, 61, 62) in terms of the DN coordinates 
one can easily find the separation relations 



H{Xi) = Xi^ + exp fii - cXi^ 



I = 



1,2, 



leading to the solution by quadratures of the Hamilton- Jacobi equations for Hi 
and H2. 

We notice that the "change of variables" (aj, 6j) t— > (Aj, /ij) is not the lift of a 
point trasformation on the configuration space; thus, there is no contradiction 
with the results of 0, stating that it is impossible to separate the 3-particle 
Toda lattice with point tranformations. 



7 Separability of Gel'fand— Zakharevich systems 



In this section we will generalize (and give proofs of) the results of the previous 
section to the case of corank k. As we will see, the picture outlined in the pre- 
vious section still holds good. The only relevant difference concerns the Stackel 
separability, which is no longer valid in general, but requires an additional as- 
sumption on the Hamiltonians. 

We consider a bi-Hamiltonian manifold (M, P, P') admitting k polynomial 
Casimir functions of the Poisson pencil P\ = P' — XP, 



such that rii + n2 + ■ ■ ■ + Uk = n, with dimM = 2n + k, and such that the 
differentials of the coefficients H^""^ are linearly independent on M. The hI"'\ 
for a fixed a, form a bi-Hamiltonian hierarchy and, in particular, Hq'^^ (resp. 
H^J) is a Casimir of P (resp. P'). We assume that the corank of P is exactly 
k, so that the Hq^\ for a = 1, . . . , k, are a maximal set of independent Casimirs 
of P. The collection of the n bi-Hamiltonian vector fields 



=PdHr = P'dH\''\, z = l,...,n„, fc = l,...,a, (7.2) 



associated with the Lenard sequences defined by the Casimirs if*-"^ is called the 
Gel'fand-Zakharevich (GZ) system^ or the axis, of the bi-Hamiltonian manifold 
M. Since standard arguments from the theory of Lenard-Magri chains show 
that all the coefficients h["'^ pairwise commute with respect to both {■, ■} and 
{■,■}', we have 




a = 1, . . . , , 



(7.1) 



i=0 
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Proposition 7.1 Let S be a symplectic leaf of P, that is, a 2n-dimensional 
suhmanifold defined by Hq^^ = ci, . . . , Hq''^ = Ck- Then the vector fields xl"'^ of 
the Lenard sequences associated with the polynomial Casimirs ( \7.J[) of {■, - j^^ on 
M define a completely integrable Hamiltonian system on S . 

We call the family {H'f'' | i = 1, . . . , Ua-, = 1, . . . , a} of the restrictions to S of 
the coefficients of the if^"-* the GZ basis of the symplectic leaf S. The lagrangian 
foliation defined by the GZ basis will be referred to as the GZ foliation of S. 

In the following subsection we will give sufficient conditions so that a sym- 
plectic leaf S of P inherits an uN structure from the bi-Hamiltonian structure of 
M. Then we will come back to the integrable system described in the previous 
proposition and we will discuss its separability in DN coordinates. 

7.1 The induced uN structure 

Our strategy to induce on a symplectic leaf S* of P a second Poisson bracket 
which is compatible with the "canonical" one is based on the geometrical con- 
siderations already mentioned at the end of Section 2. We suppose that there 
exists a A;-dimensional foliation 2 of M such that 

CI) the foliation Z is transversal to the symplectic foliation of P; 

C2) the functions that are constant on Z form a Poisson subalgebra with respect 
to both {-, ■} and {■, ■}'. 

Thus S has a (projected) bi-Hamiltonian structure. The projection of {■, ■} 
coincides with the symplectic structure {■, -j^ of S, while the projection of 
{■,■}' defines a second Poisson bracket {-j-j^ on S. Since the compatibility 
between {■, -j^ and {■, -j^ is guaranteed by the fact that the whole pencil {■, -j;)^ 
is projectable on S, we have endowed 5* with an ujN structure. We will suppose 
it to be a regular ujN manifold, in order to apply (in the open dense set where 
the eigenvalues of are distinct) the results of Section ^ and |, leaving the 
discussion of the problem of finding the conditions on (M, P, P') and Z ensuring 
the regularity of 5* for future work. 

Let (Zi, . . . , Zk) be local vector fields spanning the distribution tangent to Z. 
Because of the transversality condition, we can always normalize these vector 
fields with respect to the Casimirs Hlf' of P: 

Z,(i7('^)) = 5t . (7.3) 
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In terms of these generators, the projectabihty condition takes a very concise 
form, as shown in 

Proposition 7.2 

1. The normalized vector fields Za locally generating Z are symmetries of P, 

LzAP) = 0, (7.4) 

and satisfy 

Lz^P' = J2y>Z,, (7.5) 

b 

where = Pd{Za{Hf^)) = [Z,,P' dH^^^] = [Z,,xf\ 

2. Viceversa, suppose that there exists a k-dimensional integrable distribution 
on M which is transversal to the symplectic leaves of P and such that ( \7.4 ) 



and ( \7.^ hold for a suitable local basis {Zi, . . . , Z^) of the distribution (and for 
suitable vector fields Y^). Then the integral foliation of the distribution satisfies 
the projectability requirements CI) and C2), so that every symplectic leaf of 
P becomes an ujN manifold. Moreover, if the Za are normalized, then they 
commute. 

Proof. First of all, we recall ([0, p. 54) that the condition that the functions 
constant along Z form a Poisson subalgebra with respect to {■, ■} is equivalent 
to the assertion that the following equations hold, 

k 

Lz.P = J2<^Z^^ (7-6) 

6=1 

for some vector fields W^. This entails the validity of assertion 2, except the 
commutativity of the vector fields Za, normalized according to (|7.3|) , that can 
be proved as follows. The integrability of the distribution implies that there are 
functions 0^^, such that 

k 

[Za,Z,] = J2<Pl,Z,, 
c=l 

and evaluating this relation on the Casimirs Hq'^^ of P we easily see that (p'^j^ = 0. 

In order to prove assertion 1, we notice that the vector fields Wa are not 
uniquely defined, and can be taken to be tangent to the symplectic leaves of P. 
This is accomplished by changing 

k 
c=l 
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Indeed, 

k k k k 

6=1 c=l 6=1 6,c=l 

k 



6=1 



since Lz^idH^ , P (IHq ) = and (|7.6| ) implies that 

Thus the vector fields Wj^ in (|7.6| ) can be chosen in such a way that WI^{Hq^^) = 

(c) 

0. Now, deriving the relation PdH^ = along Za one obtains that the normal- 
ized vector fields W^^ vanish, so that, indeed, the vector fields Za are symmetries 
of P. 

As far as the second Poisson tensor P' is concerned, in the same way we can 
show that there exist vector fields tangent to the symplectic leaves of P such 
that 

k 

Lz^P' = Y,Y^AZ,. (7.7) 

6=1 

By deriving the relation P'dHlf* = x[^^ with respect to Za, one has that 

Ya = [Za,xi^^] = LzSPdH^^) = Pd{Za{Ht^)) . 
This completes the proof. 

QED 



In the sequel we will always suppose that the normalization conditions ( |7.3| ) 
on the transversal vector fields Za and the tangency conditions on the Ya are 
satisfied. For the sake of simplicity, we will also assume that the Za are defined 
on the whole manifold M, or at least in a tubular neighborhood of S. Next we 
give a useful formula for the (second) reduced Poisson bracket on 5*. 

Proposition 7.3 Let f , g be functions on a symplectic leaf S of P, and F, G 

arbitrary extensions of f , g to M. Then 

{f,g}s = {F,G} (7.8) 

k 

{f,9ys = {F,Gy + Y,{xl^\nZaiG) - xl^\G)ZaiF)) , (7.9) 

a=l 

where xf"^ = PMi^i"^ = ■}'. 
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Proof. The symplectic leaf 5* is given by the equations Hq = c", for a = 
1, . . . ,k, where the are suitable constants. The first formula simply says that 
{■, -js corresponds to the symplectic structure of S. The second formula follows 
from the remark that F := F — X]a=i ■^a.{F) (^H^"'^ — coincides with F and 
fulfills Zb{F) = on S*. Hence it can be used to compute {/, S'j's, giving { \l.9\) . 

QED 



Remark 7.4 The projectability conditions we have imposed in order to endow 
a fixed symplectic leaf 5* with an ujN structure can be weakened in the following 
way. We can consider a distribution transversal to TS and defined only at the 
points of S, generated by a family of vector fields {Zi, . . . , Z/^), normalized 
as Za{HQ^) = {dHQ\Za) = S^. Then we introduce, according to ( [7.9| ), a 
composition law {-, -j'^ on C°°{S) and we look for conditions ensuring that it is 
a Poisson bracket, compatible with {•, One can show |]T^ that {■, - j'^ is a 
Poisson bracket if and only if 



E xf'^) A (izSn + Y.[Za. XI] ^z^+\Y, X^^^ A X? A z, 

a=l \ 6=1 / a,6=l 







(7.10) 

at the points of S. In this case, the two Poisson brackets are compatible if and 
only if 

k 

Y,X['^ MzSP) = ^ (7.11) 

a=l 

at the points of S. Hence, the requirements (|7.4|) and ( [r.5|) , on the whole 



manifold M, are very "strong" solutions for (|7.1CI| ) and ( 7.11 ). Finally, we 



mention that the reduction process presented in this remark does not fit in the 



Marsden-Ratiu scheme [30|, whereas the one based on CI) and C2) clearly does. 



7.2 Separability and the control matrix 

After endowing any symplectic leaf S of P with an uN structure, we can recon- 
sider the GZ foliation of S and prove its separability in DN coordinates. Notice 
that (see also below) the restrictions to S of the bi-Hamiltonian vector fields 
X-''"'* are not bi-Hamiltonian with respect to the ujN structure of S. This is due 
to the fact that this structure is obtained by means of a projection, while the 
Hamiltonian are restricted to S. 



31 



We suppose that {Zi, . . . , Zk) are vector fields on M, fulfilling the hypotheses 
of part 2 of Proposition |7i^ and normalized, i.e., ZaiH^"^) = 6^. Then the 



expressions ( [7781) and ( |7.9|) of the reduced Poisson brackets immediately show 
that the restrictions of hI""^ to S are in bi-involution. Therefore, they are 
separable in DN coordinates. 

Theorem 7.5 The GZ foliation of S is separable in DN coordinates. 



Using once more Theorem 3J-, we can conclude that the distribution tangent 
to the GZ foliation is invariant with respect to the recursion operator A^. We 
are going to describe the form of the associated control matrix, which will be 
needed to discuss the Stackel separability of the GZ basis. 

Let g be any function on S and let G be an extension of g to M. Using ( |7.9| ) 
and the Lenard relations on the H^°'\ we have 



(7.12) 



{Ht\9ys = {Ht\Gy + J2[xi'\Ht^)Z,iG) - x[''\G)Z,{Ht^] 

b=l 

= {Hi%G}-J2z,{Hh{Hl'\G} , 

6=1 

where we have put Hna+i •= 0- Therefore, for all g G C°°(S'), 

6=1 

or, in terms of the (reduced) Poisson tensors Ps and P'g, 

P's dHl"^ = Ps dHll\ - J2 Z^^)Ps dH? . (7.14) 

6=1 

Hence, we can conclude that 

k 

N* dHl"^ = dH^l\ - Zh{Hl''^)dH['^ (7.15) 

6=1 

and read the form of the control matrix F associated with the GZ basis. Indeed, 
if we order the n functions of the GZ basis as 

Hi'\Hi'\...,HS,Hr\...,Hi'J , (7.16) 
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then we reahze that F has & k x k block form, 



JF I Ci^2 ■ ■ ■ Cl,k 

C2,l ^2 ■ ■ ■ C2,k 

Ck,i 

with Ta an na x Ha square matrix of Frobenius type of the form 



(7.17) 



•J a 



(a)^ 



1 







Za{Ht^) 1 



1 





(7.18) 



and Ca^b a rectangular matrix with Ua rows and rih columns where only the first 
column is nonzero: 



-Zb{H[ 











Zb{H^^^) ... 



z,{hI^}) ... 



(7.19) 



Remark 7.6 The vector field Y , defined in Section ^ as the Hamiltonian vector 
field associated with — |tr by the first Poisson structure, can be obtained in 
the present setting by restricting to 5* the vector field Y^^a=i ^a- Indeed, 

i^y: = pd(r£z^{Ht^) 

a=l \a=l 

and using (^) we have Y1=i Za{H^i^) = -ti F = -ifr A^. 

Thus, we have seen that GZ systems on bi-Hamiltonian manifolds admitting 
a suitable transversal foliation provide examples of non trivial (but still some- 
what special) Hamiltonian systems for which the separability condition in DN 
coordinates holds, that is, they provide interesting examples of control matrices, 
discussed in Section ^. Such matrices were introduced (in the specific example 
of a stationary reduction of the Boussinesq equation) in [|T^. 
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7.3 Stackel separability of GZ systems 

Let us now consider the Stackel (i.e., linear) separability of GZ systems. We 
have seen that the invariance with respect to N of the Lagrangian distribution 
defined by the restricted Hamiltonians H^'^'^ is a consequence of the Lenard 
recursion relations on M, and that the nontrivial coefficients in F are given by 
the deformations of the polynomial Casimirs along the normalized generators 
Za of the foliation Z. On the other hand, in Section ^ we have proved that, 
in the corank 1 case, the control matrix F automatically satisfies the condition 
for Stackel separability, N*dF = FdF. The next proposition shows that, in 
order to ensure this condition in the general case, one has to require that the 
Hamiltonians h\°''^ be affine with respect to the vector fields Za- 

Proposition 7.7 The GZ basis, formed by the h\'^\ is a Stackel basis (i.e., it 
is Stackel separable in DN coordinates) if and only if Zb{Zc{Hj'^^)) = on S, 
for all b,c,d = 1, . . . ,k and for all j = 1, . . . , n^. 



Proof. Stackel separability is equivalent to N*dF = FdF, where F is the 
control matrix ( [7.17| ). Since dF has nonvanishing entries only in the columns 
1, ni + 1, n2 + 1, . . . , Uk-i + 1, this condition takes the form 

d{z'^^)) = d{Z^l)) - J2 , (7.20) 

c=l 

where, as usual, we have put Hn^+i '■= 0. In order to compute the left-hand 
side of (|7.20|) , we observe that (|7.9| ) implies 



P'sdf= (^P' + 5^Z,AX;^)j dF , 

where / G C°°{S) and F is any extension of /. Moreover, we have that 

/ k \ k 

Lz. Ip' + Y,ZcA ) = Y,[Za, Z,] A X['^ = , (7.21) 



c=l / c=l 
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since the Z}, commute. Hence, 

k 



c=l 



= Lz(pdH\l-j2Zc{Hf^)PdHf 



c=l 



c=l 



(7.22) 



c=l 



c=l 



(7.23) 



so that 

iv* d{z'^'^)) = d{z^i)) - J2 z'^^mz'^^)) 

c=l 

-J2z^iUt^l'^))dHi^K 

c=l 

A comparison with ( |7.2CI| ) completes the proof. 

QED 

Thus, the GZ basis is Stackel separable if (and only if) the second derivatives 

of the Hamiltonians along the transversal vector fields vanish. This condition is 
automatically verified in the case of corank k = 1. This "discrepancy" between 
the generic and the rank 1 case can be understood as follows. Since, by assump- 
tion, the transversal distribution Z is integrable, the tubular neighborhood in 
which it is defined is equipped with a fibered structure, in which the fibers are 
the symplectic leaves of P. The conditions 







of equations ( [7.4|) and ( 7.21|) imply that the recursion operator (to be seen, in 
this picture, as an endomorphism of the vertical tangent bundle to the local 
fibration) is invariant along all the Za- So its eigenvalues and hence its minimal 
polynomial are invariant with respect to the Za. In the case = 1, as we have 
seen in Section ^ the coefficients of the minimal polinomial are the derivatives 
of the Casimir with respect to the (single) transversal vector field Z, but this is 
not necessarily true in the higher corank case. Notice that, whenever the second 
derivatives of the Casimirs vanish, our separated variables are "invariant" with 
respect to the Casimirs, as the one considered in ED . 
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Still under the assumptions of the above proposition, the results of Section 
§ tell us how to construct the Stackel matrix and, in principle, the separation 
relations. We also know that the entries of the Stackel matrix and of the Stackel 
vector can be used (under additional hypotheses) to explicitly find the separa- 
tion coordinates, i.e., the DN coordinates. In the next section we will exploit 
the special properties of the GZ foliation in order to determine the separation 
relations and, eventually, the DN coordinates without computing the Stackel 
matrix. 



8 Separation relations for GZ systems 

Let us consider the GZ foliation (on the symplectic leaf S) studied in Subsection 
[7.2| . The aim of this section is to write, in the Stackel separable case, the 
separation relations for the Hamiltonians of the GZ basis. To simplify the 
notations, we will not use anymore the symbol " to denote the restriction to S. 

First of all, we notice that the relevant information contained in the n x n 
control matrix F is actually encoded in the kxk polynomial matrix F(A), which 
is the Jacobian matrix of the Casimirs H^°'\\) with respect to the transversal 
(normalized) vector fields Zb, that is, the matrix 



F(A) 



Zi(iJ«(A)) ■■■ Zfc(i7(^)(A)) 



•1) 



We can translate the results about separability and Stackel separability of GZ 
systems, based on the n x n matrix equations 



N*dH 
N*dF 



FdH 
FdF 



(8.2) 
(8.3) 



into corresponding equations for the polynomial matrix F(A). To this end we 
denote by H(A) = {H^^\X) , H^'^\\) , . . . ,H^^\\)Y the A;-component vector of 
the polynomial Casimir functions, and by = {h[^\ h[^\ . 



Hi^Y and 



Fi 



Zb{H'f^) the analogs of the vector ^(A) and of the matrix F(A), con- 
structed by using the coefficients instead of the full Casimir functions 
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Lemma 8.1 The polynomial control matrix F(A) satisfies the equation 

{N* - \)dH{X) = -f{X)dH^ , (8.4) 
which is the counterpart of the matrix equation ( \8.3{ ). 
Proof. The A""'~*-coefficient of the a-th row of ( p.4|) is exactly ( [7.15 ). 

QED 

In complete analogy, we obtain the "polynomial form" of the Stackel separability 
condition (pj.3|). 

Lemma 8.2 The GZ basis is a Stackel basis iff F(A) satisfies the condition 

{N* - X)df{X) = -f{X)dfi . (8.5) 

Proof. The simplest way to prove this lemma is to expand both sides in powers 
of A. We first write (p.5|) in componentwise form as 

k 

iV*dF^(A) = ArfF^(A)- J]F^(A)rf(Fi):, 

c=l 

and then expand in powers of A, getting 

k 

c=l 



which are exactly the Stackel conditions ( [7. 201) for the GZ basis. 



QED 



The following lemma shows that the eigenvalues of can be easily obtained 
from the matrix F(A). 

Lemma 8.3 The determinant of F(A) is the characteristic polynomial of F . 
In particular, it coincides with the minimal polynomial A (A) of the recursion 
operator N , that is, 

det F(A) = det(A/ - F) = A(A) . (8.6) 

Proof. Let Aj be an eigenvalue of F. Then one can check that the relative 
(left) eigenvectors have the form 

/ 1 i \ni—2 i i \n2—l \nfc— 1 i\ 
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where := {a{, . . . , al) is a nonzero vector such that crjF(Aj) = 0. This shows 
that det F(Aj) = 0. Since det F(A) is a monic degree n polynomial and the Aj 
are distinct, we can conclude that (|8^ ) holds. 

QED 

The next step is to introduce the adjoint (or cofactor) matrix F^(A), satis- 
fying the equation 

FV(A)F(A) = F(A)F^(A) = det F(A)/ . (8.7) 

We will show that the rows of F'^(A), after a suitable normalization, provide 
Stackel function generators and play the role of the Stackel matrix. If o"(A) : = 
Cfc F'^(A) is a row of the adjoint matrix, then, obviously, 

a(A) F(A) = A(A)e, . (8.8) 

Let crj(A) be a nonvanishing entry of ct(A) and let us consider the normalized 



row 



which satisfies the equation 



p{X) = ^a, (8.9) 
(Tj(A) 



p(A) F(A) = . (8.10) 



Proposition 8.4 Suppose that the component Pa{X) of p{X) is defined for A = 
Xi, i = 1, . . . ,n. Then it is a Stackel function generator, that is, it verifies the 
equation 

(N* - X)dpaiX) = 0, forX = Xi,i = l,...,n. (8.11) 



Proof. It is convenient to consider the full vector p(A). From equation ( |8.1U| ) 
we have 

(A^* - A)rfp(A) ■ F(A) + p(A) • {N* - A)rfF(A) = for A = A^ . (8.12) 
Using Lemma |8.2| we can write the second summand in this equation as 

p{X) ■ {N* - X)dF{X) = -p(A)F(A)rfFi for A = A^, (8.13) 
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so that we finally obtain 

{N* - \i)dp{\i) -FiXi) = for z = l,...,n. (8.14) 

But the kernel of F(Aj) is 1-dimensional, due to the fact that the Aj are distinct. 
Indeed, from ( p.7| ) we have that 

n 

det F^(A) = (det F(A))*^-^ = JJ(A - X,)''-^ . (8.15) 

i=l 

If dimker F(Aj) > 2 for some i, then the rank of F(Aj) would be less than k — 1, 
so that F^(Aj) = 0, and therefore F^(A) = (A — Aj)F(A) for some polynomial 
matrix F(A). But then det F^(A) = (A - Ai)''det F(A), contradicting (g.l5|) . 
Coming back to (^.14|) , we can assert that there exist 1-forms z/j such that 

(N* - \i)dpiK) = UipiXi) for ^ = l,...,n. (8.16) 

Since the j-th component of p(Aj) is 1, we have that the Ui vanish, and this 
closes the proof. 

QED 

Now we are ready to show how to compute the separation equations for GZ 
systems. 

Proposition 8.5 Let the Pa{X) be as in the previous proposition and suppose 
that they are defined for A = Aj. Then Ylt=i '^aWH^°'\X) is a Stdckel function 
generator. 

Proof. Let us write compactly Yl'l=i piX)H^'^\X) = p(A) ■ ^^(A) and compute 

{N* - X)d{p{X) ■ H{X)) = {N* - X)dp{X) ■ H{X) + p{X){N* - X)dH{X) . 

(8.17) 



For A = Ai the first summand vanishes thanks to Proposition while the 
second equals (according to Lemma |8Tl| ) 



-p(A,) F(A,) dH, , 

and so vanishes as well. 

QED 

Therefore we have shown that the separation relations of the GZ basis (in the 
Stackel case) are given by 

k 

J2Pa{\)H^''\\) = M^^,^^^) , ^ = 1, • • • , « , (8.18) 
0=1 
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that in the corank 1 case boils down to equation 

We end this section with the following remark. Let us suppose that the mul- 
tipliers Pa(A) and the coordinates /ii, . . . be related by a "simple" algebraic 
expression, e.g., that there exist integer numbers Pi = 0, p2, ... ,Pk such that 

/if" = Pa{K) 5 i = 1, . . . ,n and a = 1, . . . , k . 

This means, according to the results of Section |^, that the Pa^th root Pa(A) is a 
Stackel function generator satisfying the equation (|5.9| ), i.e., Y{ Pa (A)) = 1, 
for A = Aj. Then the separation relations ( p.l8| ) "degenerate" to a single one, 
that is, they can be read as the vanishing of the two-variable function 

k 

^p^"if(")(A)-<l>(A,p) (8.19) 

a=l 

evaluated at the points (Aj, /ij), for z = 1, . . . , n. Hence, in such an instance, we 
can associated with the GZ system a "spectral curve" over which the separation 
coordinates lie. 

This is an indication, which is verified in several concrete examples, that the 
theory herewith presented may provide an effective bridge between the classical 
theory of the Hamilton- Jacobi equation and its modern outsprings, related to 
algebraic integrability. In this respect, several questions naturally arise, namely, 

1. Can the degeneration property of the separation relations be characterized 
in terms of the bi-Hamiltonian structure? 

2. In this case, what can one say about the algebraicity of the separation 
relation? 



We will further address these problems in [|T0[. In this paper we limit ourselves 



to give an example related to loop algebras, where all these features are present. 

9 An example related to 5i{3) 

Applications of the scheme we have presented in this paper have already ap- 
peared in the literature. Namely, in a preliminary picture of these ideas 
has been applied to the ts-stationary reduction of the Boussinesq hierarchy. 



Subsequently, in [11| we have shown how to frame all stationary reductions of 



the KdV theory inside this picture, and in |12] the classical A^-Toda lattices 



have been considered (see also PB[ for the Neumann system, and M). In this 
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final section we will illustrate how our theoretical scheme concretely works in an 
example, which is related to the ts-stationary Boussinesq system, in the sense 
that the latter can be obtained via a reduction from the one we will present. 
Even if, for the sake of brevity, we will stick to such a particular example, we 
claim that the same arguments hold for a wide class of integrable systems on 
finite-dimensional orbits of loop algebras, studied, e.g., in In these cases, 
the DN (separation) coordinates turn out to be the so-called spectral Darboux 
coordinates |[. 

The system we are going to study is defined on the space s[(3) x sl{3) of pairs 
(Xo,Xi) of 3 X 3 traceless matrices. The cotangent (and the tangent) space at 
a point is identified with the manifold itself via the pairing 



((^0, V,), {Wo, m)) = tr (Wo + V^W_ 



1) , 



so that the differential of a scalar function F is represented by a pair of matrices, 

,^ , OF OF 
dF 



dXo dXi 



We introduce |]28|, |3^ the two compatible Poisson tensors defined, at the point 
(Xo,Xi),by 



P : 





1—^ 







[Xi,K)] + [Av^i] 



where 











P' : 















A 





1 
1 



One can easily see that the functions 



Ci(Xo,Xi) = tr(AXi) = (Xi)i2 + (Xi)23 
C2(Xo,Xi) = tr(A2Xi) = (Xi)i3 
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are common Casimirs of P and P'. Thus the bi-Hamiltonian structure can be 
trivially restricted to 



M = {(Xo,Xi) G s[(3) X 5[(3) I (Xi)i2 + (Xi)23 = , {X,),^ = 1} , 
which is the 14-dimensional manifold where our GZ system will be defined. 



Indeed, it can be directly shown (see also fS^]) that, if 

L(A) = X^A + AXi + Xo , 

then 



if« = -trL(A)2 and H^^^ = -tr LiX)"^ {9.1] 
2 3 



are Casimir functions of the Poisson pencil P\ = P' — XP. One finds that 
if« = A3 + i7«A2 + /7WA + /7^') 



= A^ + H^h' + Hf^X' + Hfh' + Hfh + Hf^ 



(9.2) 



where 



H^^~^ = tr {AXo + iXi^) i/f ) = tr (XqXi) i/^'^ = |tr Xq' 
i/i') = tr {A^Xo + AXi^) Hf^ = tr ^X^^ + AXqXi + AXiXq) 
i/f ) = tr (Xi^Xo + AXo') ^ = tr(XiXo2) = |tr Xq' . 

Obviously, Hq^^ and ifo'"* are Casimirs of P, whereas H2^^ and H^'^ are Casimirs 
of P'. Since the differentials of the functions H^'^^ are linearly independent on 
a dense open subset of M, and the corank of P and P' is 2, we can conclude 
that the hypotheses of Section ^ are verified, with A; = 2, ni = 2, and n2 = 4. 
The GZ system on M is given by the 6 bi-Hamiltonian vector fields associated 
with the coefficients of the Casimirs ( p.2| ). The first vector fields of the two 
bi-Hamiltonian hierarchy are, respectively, 

X« = ([AXo], [A,X,]) , Xf) = ([A^Xo], [A^Xl]) . (9.3) 

Let us fix a symplectic leaf 5* of P, defined by the constraints Hq^^ = c^, 
Hq^^ = (? . According to Proposition |7. 1| , the 6 remaining Hamiltonians define a 
Lagrangian foliation, called the GZ foliation, on the 12- dimensional symplectic 
manifold 5*. The results of Section |^ entail that, in order to separate the GZ 
system, we need a distribution which is transversal to the symplectic leaves 



42 



of p. More precisely, let {■,•} and {■,•}' be the Poisson brackets associated 
with P and P'. Then we must find a pair of vector fields {Zi, Z2), spanning a 
2- dimensional integrable distribution on M, such that 



5: 



(9.4) 



and such that the functions invariant along the distribution form a Poisson 
subalgebra with respect to both {■, ■} and {■,■}'• It is not difficult to show that 
these requirements are fulfilled by 



Z2 



E23 
Ei3 



Xi = 
Xi = 



where Eij is the matrix with 1 in the (i, j) entry and elsewhere. In fact, a 
function F G C°°(M) is invariant with respect to both Zi and Z2 if and only if 

and such functions form a Poisson subalgebra, because the matrices fulfilling 
( |9.5| ) are a Lie subalgebra of 51(3). Moreover, Zi{Hq^^) = ti {AE23) = 1 and 



tr (y4^£'23) = 0, and the analogous equations for Z2 hold, so that we 



have the normalization (|9.4| ). Then Proposition implies that 

2 

Lz.P = 0, Lz^P' ^ 



6=1 



with 



[Za,x['''^]. The vector field (see Remark |7.6|) 

Y = Y,' + Yi = [Z,,x['^] + [Z2,x[') 
is given, on account of (|9.3|), by 



Xn 



-10 
0-10 
2 



Xi = 



Hence, the symplectic leaf S has an uN structure and Theorem tells us that 
the above-defined GZ foliation is separable in DN coordinates. 

Now we will use the results of Section ^ and to discuss the Stackel separa- 
bility and the separation relations of the GZ basis. Indeed, H^^\X) and if*^^-*(A) 
are easily seen to be affine with respect to the transversal vector fields. 



ZaiZ,{H^'\\))) = Za{Z,{H^^\\))) = for all a, 6 = 1,2, 
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meaning that the GZ basis is Stackel separable. 

A set of special DN coordinates {Xi, f^i)i=i,...fi on S is determined as follows. 
We write the polynomial matrix 



F(A) 



Zi(i7W(A)) Z2(i/«(A)) 



32 



L(A) 



31 



(L(A)2)32 (L(A)2)3i 



whose determinant gives the minimal polynomial of the recursion operator 
of S: 



det F(A) = A*^ - ^ aX 



6-i 



(9.6) 



i=l 



Its roots are the eigenvalues (Ai, . . . , Ae) of A^. The /ij coordinates can be found 
with the strategy described in Section |, which consists in looking for a Stackel 
function generator /(A) such that Y{f{X)) = 1. In Section ^ we saw that the 
normalized rows of the adjoint matrix F^(A) of F(A) provide Stackel function 
generators. We have 



F-(A) 



(L(A)2)3i -L(A)3i 

-(L(A)2)32 L(A)32 



so that /(A) := — (L(A)^)3i/L(A)3i is a Stackel function generator. Since 

F((L(A)2)3i) = -L{X)si and Y{L{Xhi) = , 
we obtain Y{f{X)) = 1, and therefore 



= fix,) = -(L(A,)')3l/^(Ai) 



31 ; 



1,...,6 , 



(9.7) 



form with the Aj a set of special DN coordinates. 

At this point we could, in principle, use ( |9.6| ) and ( |9.7| ) to explicitly write 
the point {Xq,Xi) of S in terms of (Aj, /ij)j=i^...^6; and we could compute the 
functions $j in ( ^.18|) in order to obtain the separation relations for the GZ 
basis: 

Pi(A,)if«(A,) +p2(A.)/f(')(A.) = , 
with pi(A) = /(A) and P2{X) = 1. Thus we have 

fx,H^'\X,) + H^^\X,) = $.(A„/i,) . (9.8) 



44 



However, we can directly show that these separation relations coincide with the 
ones given by the spectral curves, i.e.. 



det(/i/-L(A)) = 



Since det(/i/ — L{X)) = /i^ — |tr {L{X)^)jj, — |tr (L(A)^), the points (Aj, /ii)j=iv,6 
given by (|9.6|) and ( p.7|) belong to the spectral curve if and only if 



for all A, such that 



(9.9) 



)31 



L(A,)3i(i:(A,)2)32 = 0. 



(9.10) 



Since it can be checked that equation (9^) holds for every traceless 3x3 matrix 
fulfilling (|9.10 ), we have indeed shown that the separation relations ( |9.8|) are 
given by the spectral curve. 
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